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ABSTRACT 

We studied the properties of the intra-cluster medium (ICM) in two clusters of galaxies 
(AWM 7 and Abell 1060) and two groups (HCG 62 and NGC 507) with the X-ray Observatory 
Suzaku. Based on spatially resolved energy spectra, we measured for the first time precise cu- 
mulative ICM metal masses within 0.1 and ~ 0.3 rigo. Comparing our results with supernova 
nucleosynthesis models, the number ratio of type II (SNe II) to type la (SNe la) is estimated to 
be <~ 3.5, assuming the metal mass in the ICM is represented by the sum of products synthesized 
in SNe la and SNe II. Normalized by the K-band luminosities of present galaxies, and including 
the metals in stars, the integrated number of past SNe II explosions is estimated to be close to or 
somewhat higher than the star formation rate determined from Hubble Deep Field observations. 

Subject headings: clusters: general — clusters: individual(Abell 1060, AWM 7, HCG 62, NGC 507) 



1. Introduction 

The elemental abundance of the intra-cluster 
medium (ICM) gives us important clues for under- 
standing the chemical history and evolution pro- 
cess of clusters of galaxies. Because clusters are 
the largest gravitationally bound systems in the 
universe, they are expected to confine all the met- 
als provided by member galaxies since the cluster 
formed. The amount of each metal in the ICM is 
almost equal to the integrated sum of that syn- 
thesized by supernovae type la (SNe la) and type 
II (SNe II). The expected elemental products of 
each type have been intensively studied by several 
authors (Iwamoto et al. 1999; Nomoto et al. 2006). 
X-ray measurements of the mass of each element in 
the ICM enable us to examine the past occurrence 
numbers and the ratio of SNe II to SNe la. These 
numbers also reflect the past stellar initial mass 
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function (IMF) and star formation rate (SFR) in 
clusters. 

ASCA first measured metal distributions in the 
ICM (Mushotzky et al. 1996; Fukazawa et al. 1998, 
2000), and Dupke & White (2000) and Baum- 
gartner et al. (2005) investigated the contributions 
from SNe la and SNe II to the ICM. Recent X- 
ray observations with XMM-Newton determined 
the spatial distribution and elemental abundance 
pattern of the ICM metals (Matsushita et al. 2003, 
2007b; Tamura et al. 2004; Bohringer et al. 2005; 
de Plaa et al. 2006, 2007; Werner et al. 2006). 
However, the abundance profiles in the outer re- 
gions of clusters are still poorly constrained, espe- 
cially for O and Mg which are major SNe II prod- 
ucts, because of the intense intrinsic Al-K back- 
grounds and the non-Gaussian line profiles of pre- 
vious missions in the energy range below 1 keV. 
Based not only on the low and stable background 
but also on the good sensitivity to emission lines 
below ~ 1 keV of the Suzaku XIS (Koyama et al. 
2007) , a reliable determination of O and Mg abun- 
dances to cluster outer regions has become feasible 
(Matsushita et al. 2007a; Sato et al. 2007a). 

We use H a = 70 km s" 1 Mpc" 1 = 0.7 h wo , 
fl\ = 1 — Qm = 0.73, and the virial radius 
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is assumed ri 80 = 1.95 h^ Q y/k(T) /10 keV Mpc 
(Markevitch et al. 1998) in this paper. Errors are 
90% confidence region for a single interesting pa- 
rameter. 

2. Metal mass determination 

We selected two clusters, AWM 7 and Abell 1060 
(hereafter A 1060), and two groups, HCG 62 and 
NGC 507, that were observed with Suzaku. De- 
tailed information and the analysis method for 
each object are described in Sato et al. (2007a, b,c) 
and Tokoi et al. (2007). We excluded Ni and Nc 
abundances in our analysis in this paper since 
these two elements were not reliably determined 
due to the strong and complex Fe-L line emis- 
sions. For the first time, we determined the radial 
abundance profiles of O, Mg, Si, S and Fe, out 
to ~ 0.3 ri 80 for each object. Note that the O 
abundance in the ICM is strongly affected by the 
foreground Galactic emission, because we can- 
not resolve the cluster and local lines. Combin- 
ing the abundance profile obtained with Suzaku 
and the X-ray luminous gas mass profile with 
XMM-Newton, we calculated cumulative metal 
mass within 0.1 ri 80 and the whole observed re- 
gion (table 1). Errors of the metal mass plotted 
in figure 1 are taken from the statistical errors of 
each elemental abundance in the spectral fits at 
~ 0.1 n 8 o- 

3. Deconvolution of elemental mass pat- 
terns with SNe yields 

3.1. Mass products by SNe la and II 

In order to examine the SNe la and SNe II con- 
tribution to the ICM metals, the elemental mass 
pattern of O, Mg, Si, S and Fe was examined 
within 0.1 ri 80 and the whole observed region. 
The mass patterns were fit by a combination of 
average SNe la and SNe II yields per supernova, as 
shown in figure 1. The fit parameters were chosen 
to be the integrated number of SNe la (Ai a ) and 
the number ratio of SNe II to SNe la (N u /N u ), 
because A^ a could be well constrained due to rela- 
tively small errors in the Fe abundance. The SNe 
la and II yields were taken from Iwamoto et al. 
(1999) and Nomoto et al. (2006), respectively. We 
assumed a Salpeter IMF for stellar masses from 
10 to 50 M Q with the progenitor metallicity of 



Z = 0.02 for SNe II, and W7, WDD1 or WDD2 
models for SNe la. Table 1 and figure 2 summa- 
rizes the fit results. 

All the objects exhibit similar features. The 
abundance patterns were better represented by the 
W7 SNe la yield model rather than WDD1. The 
number ratio of SNe II to SNe la with W7 is ~ 3.5, 
while the ratio with WDD1 is - 2.5. The WDD2 
model gave quite similar results as W7. Almost 
3/4 of the Fe and ~ 1/4 of the Si is synthesized by 
SNe la, in the W7 model, as demonstrated in lower 
panels of figure 1. We also examined the SNe II 
yields with Z = 0.004 in Nomoto et al. (2006), and 
the results tended to show smaller Nu/Ni a ratios 
by ~ 30% and larger x 2 values. Somewhat earlier 
estimates of SNe II yields summarized in Iwamoto 
et al. (1999) are between Z = 0.02 and Z = 0.004, 
but our x 2 values are almost the same with those 
of Z = 0.02. The values in Table 1 imply that 
the fit to the data and the Nu/Ni a ratio behave 
in a similar manner for different supernova models 
between r < 0.1 ri 80 and r > 0.1 ri 80 . 

We note that most of the fits were not formally 
acceptable based on the x 2 values in table 1. The 
90% and 99% confidence values for 3 degrees of 
freedom arc x 2 = 6.25 and 11.34, respectively. 
The models adapted here are probably too sim- 
plified, such as assumptions of uniform IMF and 
SNe yields at all times, spherical symmetry in the 
ICM, and rather crude treatment of the Galactic 
emission, etc. For example, a flatter slope of the 
IMF than Salpeter tends to yield a larger amount 
of O. It is notable that the observed amount of Mg 
is about twice as much as the model prediction 
for AWM7 and NGC 507, similar to the feature 
reported by Werner et al. (2006b) for M 87. 

3.2. Comparison with other studies 

dc Plaa et al. (2007) also found that the num- 
ber ratio of core-collapse SNe (SNe Il+Ib+Ic) to 
SNe la (using W7) is - 3.5, based on XMM- 
Newton observations of 22 clusters within 0.2 r 50 o 
for elements between Si and Ni. They also exam- 
ined the fit including the average O and Ne abun- 
dances for two clusters observed with the XMM 
Reflection Grating Spectrometers (RGS), however 
RGS observations are limited to the central clus- 
ter region where the abundance is strongly affected 
by the cD galaxy. The number ratio of SNe II to 
SNe la they have obtained is consistent with our 
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Table 1: Integrated number of SNe I (N ia ) and number ratio of SNe II to SNe la (N u /Ni a ). 



Obiect 


z 1 k(T) 1 risn 


Region SNe la Model N u 


1 1 v la 


y 2 /dof 


AWM 7 


0.01724 


< 0.1 r 180 


WJ'1 

W 7 


1.2 ± v. 2 X ID 


4.0 ± 1.2 


12.6/3 




3.5 keV 


< 0.35 ri 80 


WJ1 

W7 


7.0 ± l.U X ID 


3.7 ± 1.2 


15.9/3 




1 65 Mdc 


< 0.35 riso 


VVJJJJ1 


y.u ± l.z X 1U 


2.3 ± 1.1 


35.5/3 






<. U.OO / 180 


VV i J \ J 


7 q _i_ i n v 1 O 9 
I.Oxl.U a 1U 


3.5 ± 1.3 


20.5/3 


A 1060 


0.0114 


< 0.1 riso 


W7 


5.1 ± 1.1 x 10 s 


3.0 ± 1.6 


4.5/3 




3.0 keV 


< 0.25 ri 80 


W7 


1.5 ±0.3 x 10 9 


2.6 ± 1.5 


8.5/3 




1.53 Mpc 


< 0.25 ri 80 


WDD1 


1.6 ±0.4 x 10 9 


2.0 ± 1.4 


16.0/3 






< 0.25 riso 


WDD2 


1.4 ±0.3 x 10 9 


2.6 ± 1.6 


10.2/3 


NGC 507 


0.01646 


< 0.1 n 80 


W7 


8.2 ± 1.7 x 10 y 


4.2 ± 1.4 


6.6/3 




1.5 keV 


< 0.24 riso 


W7 


3.1 ±0.6 x 10 8 


3.7 ± 1.2 


11.7/3 




1.08 Mpc 


< 0.24 riso 


WDD1 


3.3 ±0.7 x 10 8 


2.9 ± 1.2 


22.0/3 






< 0.24 riso 


WDD2 


3.0 ±0.6 x 10 8 


3.6 ± 1.2 


13.1/3 


HCG 62 


0.0145 


< 0.1 nso 


W7 


1.0 ±0.2 x 10 8 


3.1 ±1.3 


4.2/3 




1.5 keV 


< 0.21 n 80 


W7 


2.8 ±0.6 x 10 8 


2.6 ± 1.0 


3.0/3 




1.08 Mpc 


< 0.21 n 80 


WDD1 


2.7 ±0.6 x 10 8 


2.0 ± 1.1 


12.7/3 






< 0.21 nso 


WDD2 


2.6 ±0.5 x 10 8 


2.5 ± 1.1 


4.8/3 



results. 

Tsujimoto ct al. (1995) determined that the 
number ratio of SNe II to SNe la for our Galaxy is 
~ 6.7 based on stellar metallicity, while it is 3.3-5 
for the Large and Small Magellanic Clouds. These 
results are slightly higher than those for our clus- 
ters and groups. It is plausible that spiral galaxies 
in the Local Group maintain star forming activ- 
ity until recently. We also note that our results 
are derived for metals only in the ICM, therefore 
metals in member galaxies (including stars) which 
accounts for ~ 50% of the total are not considered. 

3.3. SNe II per unit galaxy luminosity 

Based on the luminosity function at high red- 
shift, the star formation history of the universe 
has been studied. Madau et al. (1998) derived 
a simple stellar evolution model for field galaxies 
in the Hubble Deep Field (HDF). The SFR rises 
sharply, by about an order of magnitude, from the 
present to a peak value in the range 0.12-0.17 M 
yr" 1 Mpc" 3 at z ~ 1.5. Mannucci et al. (2007) 
and Nagamine et al. (2006) also derived the cosmic 
star formation history from various observations, 
with no correction for dust extinction in cosmol- 
ogy of ft-ioo = 0.7 and f^A = 0.7. 

If we assume the maximum SFR is 0.1 M Q 



Table 2: Number of SNe II (< 0.1 n 80 ) per vol- 
ume (Nii/V), K-band galaxy luminosity per vol- 
ume (Lk/V), and Nu/Lk in HDF and clusters. 



N u /V L K /V Nu/Lk 

(Mpc' 3 ) (L Q Mpc- 3 ) (L- 1 ) 

BJW <J(F 4.0xl0 8 2.5xl0~ 3 

AWM 7* 2.2 x 10 11 5.4xl0 13 4.1xl0~ 3 

A 1060* 8.2 x lO 10 5.1xl0 13 1.6xlQ- 3 

* For AWM 7 and A 1060, the number ot SNe 11 is es- 
timated from the metals of the ICM, not including the 
member galaxies and stars. 



yr _1 Mpc^ 1 and integrate it over a Hubble time, 
the total produced stellar mass is 10 9 Mq Mpc -1 . 
Supposing a Salpeter IMF between 10-50 Mq, es- 
timated number of SNe II per unit volume is < 10 5 
Mpc~ 3 for field galaxies. When we scale our clus- 
ter results only by volume, the normalized num- 
bers of SNe II within 0.1 ri 80 for AWM 7 and 
A 1060 are 2.2 x 10 11 and 8.2 x 10 10 Mpc~ 3 , re- 
spectively. These factors of ~ 10 5 discrepancies 
are primarily caused by the difference in the num- 
ber of galaxies between the field and clusters. We 
therefore scale these results with the K-band lu- 
minosity density as summarized in table 2. Using 
an average luminosity density Lk/V = (5.74 ± 
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0.86) x 10 8 h wo L Q Mpc~ 3 in Cole et al. (2005), 
the normalized number of SNe II per unit galaxy 
luminosity is 2.5 x 10 ~ 3 Lq 1 for field galaxies. The 
galaxy luminosity densities of AWM 7 and A 1060 
are taken from the Two Micron All Sky Survey. 1 
As shown in table 2, the numbers of SNe II per 
unit luminosity in the clusters are quite consis- 
tent with the above value independently estimated 
from the star formation history. However, this es- 
timate considers only the ICM metals, and inclu- 
sion of stellar metals would roughly double the 
required SNe II rate. Even though systematic un- 
certainty in this simple estimation may account 
for most of the factor of two difference, it is inter- 
esting that X-ray metals suggest incompleteness 
in the optical detection of distant star formation. 
Also, an environmental difference between fields 
and clusters may be of relevance. 
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Fig. 1.— Fit results of each elemental mass for (a) AWM 7, (b) A 1060, (c) NGC 507, and (d) HCG 62. 
Top panels show the mass within the whole observed region (black) and within 0.1 riso (red) fit by 
[Aia {(SNe la yield) + (An/Ai a ) (SNe II yield)}]. Blue dashed and solid lines corresponds to the contri- 
butions of SNe la (W7) and SNe II within 0.1 rign, respectively. Ne (atomic number = 10) is excluded 
in the fit. Mid and lower panels indicate ratios of data points to the best-fit, and fractions of the SNe la 
contribution to total mass in the best-fit model for each element, respectively. 
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Fig. 2. — (a) Integrated numbers of SNe la (W7) and SNe II plotted against gas mass within the calculated 
region, (b) Number ratio of SNe II to SNe la with W7 model. 
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